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Abstract 

We consider the eigenvalue problem of certain kind of non-compact linear operators 
given as the sum of a multiplication and a kernel operator. A degenerate kernel 
method is used to approximate isolated eigenvalues. It is shown that entries of the 
corresponding matrix of this method can be evaluated exactly. The convergence 
of the method is proved; it is proved that the convergence rate is 0{h). By some 
numerical examples, we confirm the results. 
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1 Introduction 

Consider tlie following non-compact operator defined on X = £°°[a, 6] by 

rh 

Af{x) = a k{x,u) f{u) du — x'^ f{x), VxG[a, &]. (1) 

J a 

Here, [a,b] is a compact interval, /c(.,.) G C°([a, b] x [a, 6]), and a is a real 
constant. The non-compactness of the operator A is due to the involved mul- 
tiplication operator with the corresponding function x'^ [S]. 
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We are concerned with the numerical approximation of isolated eigenvalues of 
the following eigenvalue problem: 

Xf = Af. (2) 

Eigenproblem ([2]) has some applications in electromagnetism [5]. Also, a more 
general form problem ([2]) describes the COA model in population genetics [2]. 
In [7], this model has been treated numerically by Nystrom and Galerkin meth- 
ods in the presence of some extra assumptions. A degenerate kernel method 
with piecewise linear interpolation with respect to the variable u has been 
developed for problem ([2]) in [S]. In this method, isolated eigenvalues of prob- 
lem ([2]) are approximated by eigenvalues of an operator An with rank n; but 
the entries of the matrix eigenvalue problem associated with An can not be 
evaluated exactly. Hence, we need to replace various integrals by numerical 
quadratures. This perturbation in entries causes unacceptable approximations 
for ill-conditioned eigenvalues. 

In in a degenerate kernel method is proposed for the numerical solutions of 
integral equations of the second kind with a smooth kernel which does not 
need any numerical quadrature, and the entries of its corresponding matrix 
are evaluated exactly if we ignore round-off errors. Genaneshwar [3] has been 
extended this method to the corresponding eigenvalue problem with the same 
benefit. In this paper we extend this method for the numerical solutions of 
eigenproblem ([2]). 

The rest of this paper is organized as follows: In Section [21 we describe the 
degenerate kernel method for numerical solutions of problem ([2]). Convergence 
of the method and its rate are discussed in Section [2l In Section HI we give 
some numerical results to illustrate the accuracy of the method. 

We recall that X = C°°[a, b] throughout this paper. 



2 Degenerate kernel method 

We begin this section by some preliminaries. 

Let Y be any Banach space over the complex field C. We denote the space 
of bounded linear operators from Y into Y by BL(Y). Let T G BL(Y). The 
resolvent set of T is given by 

p{T) = {zeC : (T - ziy^ e BL(Y)}. 

The spectrum of T, denoted by cr{T), is defined as cr(T) = C\p(T). The point 
spectrum of T consists of all A e cr(T) such that T — XI is not one-to-one. 
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In this case A is called an eigenvalue of T. If A is an eigenvalue of T, then 
the smallest positive integer I such that ker(T — A/)' = ker(T — XlY^^ is 
called the ascent of A. The dimensions of ker{T — XI) and ker{T — AJ)' are 
called geometric multiplicity and algebraic multiplicity of A, respectively. If 
the algebraic multiplicity of A equals one, then it is called a simple eigenvalue 
of T. 

We now describe the degenerate kernel method. Define the operators K, and 
Al on X as follows: 

rb 

K, f{x) — a k{x,u) f{u) du, (3) 

Ja 

and 

Mf{x)=x'f{x). (4) 

Then A = IC-M. 

For an integer n > 1, consider the following partition of [a,b]: 

a — Xq < Xi < ■ ■ ■ < Xn-l < Xn — b. 

Let Ij = [xj-i, Xj], hj = Xj — Xj^i for j — 1,2, ... ,n. Also, let h — max hj 

denotes the norm of the partition. We assume that /i — > as n ^ oo. For a 
positive integer r, let Br — {ti, . . . ,Tr} he the set of r Gauss points, i.e., the 
zeros of the Legendre polynomial ^{t'^ — in the interval [—1,1]. Define 
fj : [—1,1] —>■ Ij as follows: 

Then S = U^^i fj{Br) = {^ij = fj{Ti) : t = 1, ... ,r, j ^ 1,2, ... , n} is the set 
of Nfi — nr Gauss points on [a, b]. For each i = 1, . . . , r, let 

(^\ _ {X-Ti)...{x- Ti.i){x - Tj+i) ...{X-Tr) ^^rin 
'^i\X) — , , , , , . , , X fc [ 1,1 , 

[Ti - Ti) . . . [n - Ti-i)[Ti - Ti+i) ...[n- Tr) 

be the Lagrange polynomial of degree r — 1 on [—1, 1], which satisfies •^^(tj) = 
5ij. It is notable that for r = 1, = 1 on [—1,1]. Define 



(t>pq{x) = 



0, otherwise. 



Then (fipqi^ij) = SpiSjg, for p, i = 1, . . . , r and j,q ^ 1,2, ... ,n. 

We set the following notations = ^ij, ip(j-i)r+i = (pij, ioi i = 1, . . . ,r 

and j — 1,2, ... ,n. Then S — {ti : i — 1,2, ... , N^} is the set of nr Gauss 
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points in [a,b]. Now define the degenerate kernel by 

kN^{x,u) = J2J2 k{ti,tj)tpi{x)'ipj{u), x,ue [a,b]. 
i=i j=i 

The approximated operator JCn^ is defined by 

rb 

K-N^ f{x) = a kiy^{x, u) f{u) du. (5) 

J a 

For the approximation of the muhiphcation operator M., define the piecewise 
constant functions S]s[^{x) on [a, 6] as follows: 

If X e Ji, then SAr,^(a;) = mj = -] , i = 1,2, Nh- 



Now define the approximate operators /Cat,^ and M-n^, ^ as follows: 

MnJ{x) = SN^{x) f{x). (6) 
Consider now the following approximate eigenvalue problem on X: 

AjVh fNn = ^JVh /jVh (x) , (7) 

where An^ = JCn^ ~ -^n^- The rank of the operator An^ is finite, so the 
eigenvalue problem ([71) is equivalent to a matrix eigenvalue problem as we 
show in the following: 

Eigenproblem ([71) expands into 

fNh (x) = a ^i^p^ ^j) M^) - (x) fNh (x) , (8) 

p=l j=l 

where 

Cj = / 'iljj{u)fN^{u) du. 

J a 

For each i = 1, 2, . . . , N^, we multiply both sides of ([8]) by 'ipi{x) and integrate 
over the interval [a, h] with respect to the variable x to obtain 

Nh 

XN^Ci = a^kijCj - niiCi, z = 0, l,...,n, (9) 
i=i 

where 

Nh 



p=i •^'^ 
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Therefore, eigenvalues of the approximate problem ([7]) are the eigenvalues of 
the matrix A = (aij) of order N^, defined as 

ttij = akij - ruiSij, i, j = 1,2, Nh. (10) 

It is easy to see that 

/ ijji^x) ipp(yX) dx = 0, ifij^p. 

J a 

Therefore, kij reduces to 

rb 

kij = k{ti,tj) / 'tpi{x)dx. (11) 

J a 

Since ■0j for i = 1, 2, . . . , A^^ are basis functions having small support on [a, b] 
with the simple structure, the integrals Jl^ipf{x) dx can be evaluated exactly. 
Values of these integrals depend on the length of their support. Thus, if Xi for 
i = 1,2, ... ,11 are equidistance points and the kernel k{., .) is Hermitian, then 
the matrix A is normal and as a result, each eigenvalue of A is well-conditioned 
(see i). 



3 Convergence and the rate of convergence 

Theorem 1 ^T] Let Y he a Banach space. Let also T and he opera- 

tors in BL{Y). Assume that the sequence converges in norm to T, i.e., 
||T„ — T|| ^ as n — > oo. Let A he an isolated point of a{T). For each 
positive e < dist{X, a(T)\{\}) , define 

An := {An G cT{Tn) : |A„ - A| < e}. 

Then for n large enough, A„ 7^ and if Xn G A^, the sequence {A„} converges 
to A. 

In order to show that the degenerate kernel method proposed in previous 
section is convergent, it is enough to show that the operators A and {An^} 
satisfy conditions of Theorem [1], i.e., they are in BL{X) and An^, converges 
in norm to A as 00. The linearity of these operators is obvious. In [5J, 

we have shown that the operator A is bounded. By a similar discussion it is 
seen that the operator ANy, is bounded, too. 

If we assume that k{., .) G ^"([0, h] x [a, h]), then ||/C — Kn,-^ \\ = Oijf) [3]. Also, 
from [5] we have \\M. — A^Af,J| = 0{h). Therefore, 

\\A-ANj=0{h). (12) 
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Table 1 

Computational results for sample problem 1. 



n 


|An(l)-A(l)| 


ratio{l) 


|An(2)-A(2)| 


ratio{2) 


10 


0.0763 




0.0916 




20 


0.0398 


1.92 


0.0460 


1.99 


40 


0.0202 


1.97 


0.0229 


2.01 


80 


0.0100 


2.02 


0.0113 


2.03 


160 


0.0048 


2.08 


0.0054 


2.09 



On the other hand, with the assumptions of Theorem [T], if A is a simple 
eigenvalue of A, then there exists a constant c such that [H Page 201] 

|Aiv„ - A| < c ll^jv, - ^11, for all large A^^,. (13) 

From (fT2|) and (fT3|) we conclude that the proposed degenerate kernel method 
of previous section is convergent with the rate 0{h) for simple eigenvalues. 



4 Numerical results 



In this section, we solve two samples of problem ([2]) using the degenerate kernel 
method, proposed in Section [2l In order to show the convergence, the problems 
are solved for different numbers of interpolation nodes. Two arbitrary positive 
eigenvalues of this problem, denoted by A(l) and A(2), are approximated by 
An(l) and A„(2), respectively. In both examples, we take r = 2 and compute 
the absolute error and the ratio of each eigenvalue. 



4-1 Example 1 

Consider problem ([2]) with the kernel 

A;(x, u) = exp(— (-u — xY), 

and the constant a = 1 on the interval [—2,2]. Computational results are 
shown in Table [TJ 
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Table 2 

Computational results for sample problem 2. 



n 


|An(l)-A(l)| 


ratio{l) 


|An(2)-A(2)| 


ratio(2) 


10 


0.1607 




0.0712 




20 


0.0793 


2.03 


0.0287 


2.48 


40 


0.0396 


2.00 


0.0124 


2.31 


80 


0.0196 


2.02 


0.0055 


2.25 


160 


0.0095 


2.06 


0.0025 


2.20 



4-2 Example 2 



Consider problem ([2]) with the kernel 

and the constant a = 1 on the interval [—4,4]. Computational results are 
shown in Table [2l 

In both examples, convergence of solutions is seen. The ratio of error in each 
step also confirms our theoretical discussion in Section [31 



5 Acknowledgements 

The authors wish to thank Masoud Amini for valuable discussions during 
the preparation of this paper. Also, many thanks to Alaeddin Malek for his 
supports as the head of applied mathematics department of Tarbiat Modares 
University. 



References 

[1] M. Ahues, A. Largillier, B. V. Limaye, Spectral Computations for Bounded 
Operators, Chapman k Hall/CRC, 2001. 

[2] Biirger, The Mathematical Theory of Selection, Recombination, and Mutation, 
Wiley Chichester, 2000. 

[3] N. Gnaneshwar, A degenerate kernel method for eigenvalue problems of compact 
integral operators, ADVCOM 27 (3) (2007) 339-354. 



7 



[4] H. Kaneko, P. Padila, Y. Xu, Supcrconvcrgcnce of the Iterated Degenerate Kernel 
Method, Appl. Anal. 80 (3) (2002) 331Y351. 

[5] H. Majidian, E. Babolian, A Degenerate Kernel Method for Eigenvalue 
Problems of Non-compact Operators with Applications in Electromagnetism and 
Continuum-of-alleles Models, submitted for Publication. 

[6] K. E. Morrison, Spectral Approximation of Multiplication Operators, New York 
J. Math 1 (1995) 75-96. 

[7] O. Redner, Discrete Approximation of Non-compact Operators Describing 
Continuum-of-alleles Models, Proc. Edinb. Math. Soc. 47 (2004) 449-472. 

[8] Y. Saad, Iterative Methods for Sparse Linear Systems, PWS Publishing 
Company, 1996. 

[9] M. Shahabadi, S. Atakaramians, N. Hojjat, Transmission Line Formulation 
for the Full- wave Analysis of Two-dimensional Dielectric Photonic Crystals, 
lEE. Proc. Sci. Meas. Technol. 151 (5) (2004) 327-334. 



8 



